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Geometric Controls for a Tethered Quadrotor UAV 


Taeyoung Lee 


Abstract —This paper deals with the dynamics and controls 
of a quadrotor unmanned aerial vehicle that is connected to 
a fixed point on the ground via a tether. Tethered quadrotors 
have been envisaged for long-term aerial surveillance with high¬ 
speed communications. This paper presents an intrinsic form 
of the dynamic model of a tethered quadrotor including the 
coupling between deformations of the tether and the motion of 
the quadrotor, and it constructs geometric control systems to 
asymptotically stabilize the coupled dynamics of the quadrotor 
and the tether. The proposed global formulation of dynamics 
and control also avoids complexities and singularities associated 
with local coordinates. These are illustrated by numerical 
examples. 

1. Introduction 

A quadrotor unmanned aerial vehicle (UAV) consists of 
two pairs of counter-rotating rotors and propellers, located 
at the vertices of a square frame. It is capable of vertical 
take-off and landing (VTOL), but it does not require complex 
mechanical linkages, such as swash plates or teeter hinges, 
that commonly appear in typical helicopters. Due to its sim¬ 
ple mechanical structure and higher thrust-to-weight ratio, 
it has been utilized for various applications such as mobile 
sensor network, and autonomous delivery systems. 

In particular, tethered quadrotors have been envisaged for 
aerial surveillance recently. This corresponds to a quadrotor 
UAV that is connected to a fixed ground point or a mobile 
vehicle via a tether, which can be used to for various 
purposes, such as supplying power to the quadrotor con¬ 
sistently, or high-bandwidth communication. This eliminates 
certain drawbacks of typical quadrotors that rely on onboard 
batteries and that use radio communication, namely limited 
flight endurance and communication bandwidth. Therefore, 
tethered quadrotors are particularly useful for several mission 
scenarios, such as prolonged aerial surveillance, communi¬ 
cation relay, or border protection with secure, high-definition 
video feed [1]. 

Stabilization and observer design for a tethered quadrotor 
have been studied in [2], [3], [4]. However, these results are 
commonly based on two simplifying assumptions that the 
quadrotor and the tether remain on a fixed two-dimensional 
plane, and the tether is always taut. As such, they may 
not be suitable for realistic scenarios where the quadrotor 
performs aggressive translational and rotational maneuvers, 
or the tether deforms due to the dynamic coupling with the 
quadrotor. 
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Fig. 1. Tethered quadrotor UAV: flexible tether is modeled as a serial chain 
of n links, and the conflguration manifold is Q = x SO(3). 


This paper is focused on eliminating such restrictive 
assumptions both in the dynamic model and control sys¬ 
tem design. First, we present a mathematical model of the 
tethered quadrotor in the three-dimensional space including 
deformation of the tether, which is considered as an arbitrary 
number of rigid links that are serially interconnected via 
ball joints. Next, a geometric tracking control system is 
designed such that the quadrotor can asymptotically follow a 
given desired trajectory in the three-dimensional space while 
controlling the tension along the tether, assuming that the 
tether is taut. Numerical simulations illustrate that the tension 
should be sufficiently large to prevent lateral vibrations of 
the tether, when applied to the flexible tether model. This 
motivates the development of another control system to 
stabilize both the quadrotor and the tether simultaneously, 
while incorporating the deformable dynamics of the fiexible 
tether explicitly. 

Another distinct feature is that the equations of motion and 
the control systems are developed directly on the nonlinear 
configuration manifold in a coordinate-free fashion. This 
yields remarkably compact expressions for the dynamic 
model and controllers, compared with those based on local 
coordinates that often require symbolic computational tools 
due to complexity of multibody systems. Furthermore, sin¬ 
gularities of local parameterization are completely avoided. 

In short, the main contributions of this paper are summa¬ 
rized as (i) a global dynamic model for a tethered quadrotor 
with fiexible tether, (ii) a geometric tracking control system 
for the quadrotor maneuvers in the three-dimensional space 
and the tension when the tether is taut, and (iii) a geometric 
control system for the quadrotor and the fiexible tether. The 
second item can be considered as a new contribution inde¬ 
pendently, and the dynamics and control including fiexibility 
of tether have been unprecedented. 



11. Dynamics of a Tethered Quadrotor UAV 

Consider a quadrotor UAV that is connected to a fixed 
point on the ground via a tether. Define an inertial frame 
whose origin is located at the pivot point where tether is 
attached to the ground. The third axis of the inertial frame is 
pointing downward along the direction of gravity. Define a 
body-fixed frame whose origin is located at the mass center 
of the quadrotor (see Figure 1). 

We approximate the tether by an arbitrary number, namely 
n, of rigid links that are serially connected by ball joints. 
The links are indexed from the ground to the quadrotor in 
ascending order, i.e., the first link is attached to the ground, 
and the n-th link is attached to the quadrotor. Let the mass 
and the length of the i-th link be G M, respectively, 

where it is assumed that the mass is uniformly distributed 
along each link. Throughout this paper, the subscript i is 
considered as an element of { 1 ,..., n}. 

The direction of the i-th link toward the quadrotor is 
denoted by the unit-vector G = {q G | ||g|| = 1}. 
The location of the outward end of the i-th link, namely 
Xi G is given by 

i 

= ( 1 ) 

i=i 

and therefore, the i-th link connects Xi-\ with xi assuming 
^0 = Osx 1 • The n-th link is attached to the mass center of the 
quadrotor such that the location of the quadrotor corresponds 
to X = Xn = ^ Let R G SO(3) = 

{i? G I R = /, det[i?] = 1} be the rotation matrix 
describing the attitude of the quadrotor, and it represents the 
linear transformation of the representation of a vector from 
the body-fixed frame to the inertial frame. Therefore, the 
configuration manifold of the presented tethered quadrotor 
is Q = (S^)^ X SO(3). 

The kinematics equations are given by 

qi = UiX qi, (2) 

R = Rn, (3) 

where uji G is the angular velocity of the i-th link 
represented with respect to the inertial frame. Without loss 
of generality, it is assumed that qi ' uJi = 0, i.e., uJi is 
perpendicular to qi. The standard dot product is denoted by 
X ' y = x^y for any x^y G in this paper. The vector 
U G corresponds to the angular velocity of the quadrotor 
represented with respect to the body-fixed frame, and the hat 
map • : ^ 5o(3) = {S e | = —S} is defined 

such that xy = X X y for any x^y G The inverse of the 
hat map is denoted by the vee map V : 50 (3) ^ 

The dynamic model of the quadrotor is identical to [5]. 
The mass and the inertia matrix of the quadrotor are denoted 
by m G M and J G respectively. It generates a thrust 

u given by 

^ = -fRes. (4) 

with respect to the inertial frame, where / G M is the total 
thrust magnitude and 63 = [0, 0,1]^ G It also generates 


a moment M G with respect to its body-fixed frame. The 
control input of the presented tethered quadrotor is (/, M). 

A. Euler-Lagrange equations 

We derive a global form of the equations of motion for the 
tethered quadrotor via Lagrangian mechanics. The material 
points on the i-th link are parameterized as Xi-i + Qqi for 
Ci G [0,/^], and the mass of the infinitesimal element dQ 
corresponds to Thus, the kinetic energy of the i-th 

link can be written as 

T = ^ ^^Wxi-i + CiQifdC 

= . 

The kinetic energy of the quadrotor is composed of the 
translational kinetic energy and the rotational kinetic energy, 

Rquad ~ 2^ll^’^ll 2^ * TU. 

The total kinetic energy is T = Tquad + 

we have Xi = Substituting this and rearranging, 

the total kinetic energy can be written as 

1 n ^ 

T = - MijQi ■ Qj + -n ■ jn, (5) 

where the fixed inertia terms Mij G M are defined as 
Mii= (m + \mi] if + V mplf, 

^ ’ p=i+i 

for any 1 < i < n, and the off-diagonal terms are defined 
for 1 < jf < i < n as 



Next, the gravitational potential energy of the i-th link and 
the quadrotor are given by 

Ui = -rriiges • (xi_i + 

^quad — '^ 5^63 ' ^n- 


Therefore, the total gravitational potential is U = Uquad + 
written as 


/ i-\ 


U = -ges 'V^rrii ^Ijqj + -kqi ^^mluqk 


i=i \j=i 
n 

= -ge-z 

i—\ 

where Mg. G M is defined as 


k=l 


( 6 ) 


Mg.=m^-mi-\- ^ Trip. 

p=i+l 

The Lagrangian is given hy C = T — U from <0 and 0- 
A coordinate-free form of Lagrangian mechanics on the 
two-sphere and the special orthogonal group SO(3) for 


various multibody systems has been studied in [6], [7]. The 
key idea is representing the infinitesimal variation of qi e S‘^ 
in terms of the exponential map: 


A ^ 

% = -T- 

de 


exp{e^i)qi = Ci x qi, 


(7) 


5=0 


for a vector G with ’ Qi =0. This guarantees that 
the infinitesimal variation is at the correct tangent space of 
the two-sphere, i.e., 6qi G T^.S^. Similarly, the variation of 
Ri is given by SRi = Rifji for r]i G 

By using these expressions, the equations of motion can 
be obtained from Hamilton’s principle as follows. 


jn = M, 


( 8 ) 


Muqi - qf Mijiij + MuWqiW^qi + qjMg.lige:^ = -qfliU, 

(9) 


i=i 


(see Appendix A). This can be rearranged as 


M{q) 


~qi 


'Gi{q,qy 


'-hqlu 


+ 

G2{q,q) 

= 

-hqju 

_Qn_ 


Qn{q,q)_ 


-lnqlu_ 


( 10 ) 


where M{q) G ■ 


M{q) = 


3nx3n 


, Qi{q^ q) G are defined as 


Mn/ 

-qlMi2 ■ ■ ■ 

-qlMin' 


— ^2^21 

M 22 I 

-qjMzn 

, (11) 

—qnMjii 

— qnMji2 

MnnI _ 


Ui \\qi\y qi + qiMgJigez. 


(12) 


Alternatively, the equation © can be rewritten in terms of 
the angular velocity as 

n 

MiiUJi - qi'^Mij{qjdij + \\LOj\fqj) - qiMgJigez = kqiU. 

i=i 

(13) 

Together with the kinematics equations 0 and ([^, these 
describe the dynamics of the tethered quadrotor. 

III. Control System Design for Taut Tether 

Designing a control system for the tethered quadrotor 
described by (§ and ( p^ is challenging as it is highly 
underactuated: there are 2n+3 degrees of freedom but only 4 
independent control inputs. In this section, we first design a 
control system for the special case when there is a single link, 
i.e., n = 1, assuming that the cable is always taut. Instead, 
the tension along the tether is also controlled such that the 
tether is stretched even if the deformation of the tether is 
included in the dynamic model. The deformation of tether 
will be incorporated later at Section IV. Throughout this 
section, the subscript 1 is removed for brevity, i.e., q = qi 
and uj = uji. 


A. Problem Formulation 

When n = 1, the equation of motion ( p3] ) reduces to 

cj - aqe^ = fdqu, (14) 


where the constants f3 are given by 

^ Mgji ^ ^ 1 

^ Mil ^ (m-h^)/i^’ Mil (m + mi/3)/i’ 

Next, we find the expression of the tension along the tether. 
Since the location of the quadrotor is given hy x = Iq, using 
^ and its acceleration can be written as 

x = lq = l{-q(:b - \\ujfq) = l{-afez - pq^u - ||wf g). 

Let A G be the internal force exerted by the tether 
on the quadrotor. Considering the free-body diagram of the 
quadrotor excluding the link, from Newton’s second law, 
mx = u F mges + A. The tension along the link, namely 
T G M corresponds to the component of the negative internal 
force —A along the direction of the link q, i.e., T = • q. 

Note that it is defined such that a positive tension T implies 
that the tethered is being stretched. By combining the above 
two equations, 

T = q - {u-\- mges — mx) = q • {u ges) + 

In short, the equations of motion for the tethered quadrotor 
and the tension when n = 1 are given by ([8]), and 

(j — aqe^ = pqu^^ ( 15 ) 

T = mgq ' esml\\uj\\‘^q ' u\ ( 16 ) 

where G denote the component of u that is 

perpendicular to u, and the other component that is parallel 
to u, respectively, given by 

u-^ = (73x3 - qq^)u = -fu, (17) 

U^\ =qq'^u= {I 3 x 3 +q'^)u. (18) 

A tracking control problem for the tethered quadrotor 
is formulated as follows. Suppose that a smooth desired 
trajectory of the direction of the link, namely qd{t) : M ^ S^, 
and the desired tension Td{t) : M ^ M are given. The desired 
direction satisfies 


qd{t) = ujd{t) X qd{t), (19) 

for the corresponding desired angular velocity ujd{t) G 
satisfying ujd{t) • qd{t) = 0. We wish to design the control 
input of the quadrotor {f^M) such that this desired trajectory 
becomes an asymptotically stable equilibrium of the con¬ 
trolled system. 

B. Simplified Dynamic Model (n = 1) 

The presented quadrotor is underactuated since the total 
thrust is always parallel to its third body-fixed axis. This can 
be directly observed from the expression of the total thrust 
given by IX = —fRe^. The magnitude / of the total thrust 
and the total control moment M are arbitrary. To overcome 
this, we first consider a simplified dynamic model where the 
















quadrotor may generate the total thrust along any direction. 
This is equivalent to designing a desired total thrust u based 
on ([Tg and without considering its attitude dynamics 
This is possible as the attitude dynamics does not directly 
appear in the dynamics of the link at ([T^. The effects of the 
attitude dynamics will be incorporated later. 

The equations of motion for the link and the tension 
given by © and ( p^ have the following structure: the link 
dynamics is controlled by the perpendicular component of 
the control input and the tension is controlled by the 
parallel component of the control input uK Therefore, u-^ is 
designed such that the link asymptotically follows its desired 
direction, and is designed for the desired tension. The 
resulting complete control input is obtained by combining 
them together. 

First, we design the parallel component. Since the tension 
is an algebraic function of it 11 at (16), it is designed as 


mII = {Td - mgq ■ 63 - ml\\uj\\'^)q, 


( 20 ) 


such that the resulting tension is identical to always. 

Next, we design for the link dynamics ^ and (15) 
such that q as t ^ oo. Control systems for the 

unit-vectors on the two-sphere have been studied in [ 8 ], 
[9]. In this paper, we adopt the control system developed 
in terms of the angular velocity in [9]. Define the tracking 
error variables, namely G as 


e, = X g, e^=u) + fud- 


the attitude of each quadrotors is controlled such that the 
direction of its third body-fixed axis. Res becomes parallel 
with —u given at ( [ 2 ^ . 

The corresponding attitude controller is similar with [5], 
[10]. The desired direction of the third body-fixed axis is 


b3. = - 



(24) 


There is an additional one-dimensional degree of freedom in 
the desired attitude, corresponding to rotation about 63 ^. A 
desired direction of the first body-fixed axis, bi^ (t) G is 
introduced to resolve it. The resulting desired attitude is 


Rc 


WhMS 


(25) 


and the desired angular velocity is obtained by Qc = 
{R^RcY G Define the error variables for the attitude 
dynamics as 

eR = ^{Rc ^ RcV, ea = Q - R^RcO,c. 

The thrust magnitude and the moment vector of quadrotors 
are chosen as 


f = -u- Re 3 , (26) 

M = en + Qx J12 

e 

- J{tlR^RSlc - R^Rc^c), (27) 


For positive constants kq^k^j G M, the normal component of 
the control input is chosen as 

M-*- = - {q ■ ujd)q - fuJd - age^}. 

( 21 ) 


Note that the expression of is perpendicular to q by 
definition. Substituting into ( p~5) ), and rearranging it with 
the facts that the matrix — g| corresponds to the orthogonal 
projection to the plane normal to qi and qf = —qi, 

U! = -kgCg - - (q ■ U!d)q - q^tOd- ( 22 ) 


In short, the control input is given by 


u = uY 


(23) 


from ( 20 ) and ( 21 ), for the simplified dynamic model. 

Proposition 1: Consider the simplified dynamic model 
described by and ( p^ , where n = 1 and u 

can be arbitrarily selected. The control input is designed 
as ( [2^ . Then, the zero equilibrium of the tracking error, 
Yqi^uj) = (O 5 O) is exponentially stable, and the tension is 
identical to its desired value, i.e, T{t) = Td{t) for any t. 

Proof: See Appendix B. ■ 


C. Full Dynamic Model (n = 1) 

The above control system for a simplified dynamics model 
is generalized to the full dynamics model that includes the 
attitude dynamics of the quadrotor. The control 

force of the full dynamic model is given by —fRes- Here, 


where e, kR^ kQ are positive constants [5]. 

Proposition 2: Consider the full dynamics model defined 
by (|^, ([^, ([^, ( p~5] ), and ( p^ , where n = 1. Control 
inputs (f^M) are designed as ^2^ and (27), where the 
desired control force u is given by 


|. Then, there exist 
e* > 0 , such that for all e < e*, the zero equilibrium of the 
tracking errors (e^, eo) is exponentially stable, and 

\imt^ooT{t) = Td{t). 

Proof: See Appendix C. ■ 


D. Numerical Examples 

Compared with the prior results where the motion of the 
quadrotor and the tether is restricted to a two-dimensional 
plane, the proposed control system guarantees that the 
quadrotor and the tether asymptotically follows a desired 
trajectory in the three-dimensional space, while controlling 
the tension along the tether. These are illustrated by a 
numerical example as follows. 

Properties of the quadrotor are chosen as m = 0.755 kg 
and J = diag[0.0043, 0.0043, 0.0103] kgm^. The total mass 
and the total length of the tether are 0.3 kg and 5 m, 
respectively. Initially, the tether is aligned along a horizontal 
direction with zero angular velocity, i.e., g(0) = ei,ci;(0) = 
0. The desired trajectory is chosen such that the quadrotor 
follows a figure-eight curve on the sphere, 

= [cos ^ COS 0, sin^, — cos^sin^]^, 

0{f) = ^ sin0.27rf, (j){f) = ^ sin0.47rt + 
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(c) Tether direction error eg (d) Quadrotor attitude error cr 

Fig. 2. Numerical results for n = 1 (animation available at http:// 
fdel.seas.gwu.edu/CDC15_Fig2.mov) 

and the desired tension is Td = 5N. The corresponding 
simulation results are presented at Figure]^ where it is shown 
that the tracking errors converge to zero. 

Next, we apply the presented control system developed 
for n = 1 into the dynamic model of a flexible tether with 
n = 5. This is justifled with the assumption that the tether 
remains taut even for the flexible tether model, if the tension 
is sufficiently large [2], [3], [4]. When computing the control 
input, the direction of the taut tether is approximated by the 
direction from the origin to the quadrotor, i.e., q = 

Numerical results when = 10 N are illustrated at Figure 
where the snapshots of the controller maneuvers, and 
the positions of the first, the third, and the last links are 
presented. While the position of the last link that corresponds 
to the position of the quadrotor follows the desired position 
relatively well, there are nontrivial lateral vibrations at the 
first link and the third link. To reduce the vibrations, the 
desired tension is increased to Td = 20 N at Figure 
However, there still exist persistent vibrations as illustrated 
by Figure 4.(b) and animation. 

IV. Control System Design eor Flexible Tether 

The control system designed at the previous section can 
excite the lateral vibration of the tether when applied to 
the flexible tether model, and it may require increasing the 
tension of the tether unnecessarily large to avoid vibrations 
in certain cases. Motivated by these, in this section, we 
design another control system for the tethered quadrotor 
while explicitly incorporating the dynamics of a flexible 
tether. For simplicity, the desired configuration is selected 
as qi^{t) = — 03 , i.e., all of the links are aligned along the 



Xl,X3,X5 = X 

Fig. 3. Numerical results for n = 5, = 10 N (animation available at 

http://fdcl.seas.gwu.edu/CDC15_Fig3.mov) 


(a) Snapshots of maneuver (b) Position of selected links 

Xl,X3,X5 = X 

Fig. 4. Numerical results for n = 5, = 20 N (animation available at 

http://fdcl.seas.gwu.edu/CDC15_Fig4.mov) 

direction of the gravity, and the quadrotor is located directly 
over the pivot point at Xd = 

A. Simplified Dynamic Model (n > 1) 

Similar to the prior section, we first consider the simplified 
dynamic model where the total thrust u can be arbitrarily 
selected. The proposed control system is composed of two 
parts: an output tracking controller to translate the quadrotor 
position X into the vicinity of Xd, and a control system to 
asymptotically stabilize the desired configuration. 

1) Tracking Control for Quadrotor Position: Here, we 
design a control input u such that the position of the 
quadrotor, namely x G is translated into an intermediate 
point {1 — S)xd for a. constant 0 < 6 < 1. Note that when 6 
is small, the intermediate point becomes closer to the actual 
desired position Xd. Also, (1 — S)xd belongs to the set, 

n 

D, = {xeRy\\x\\<Y,h}, (28) 

i=l 

which is the sphere centered at the origin whose length is 
strictly less than the total length of tether. 

From ( p^ , qi can be written as 

n 

Qi = M/ (q) {-Qj{q, q) - ljq]u) , (29) 

i=i 



















































where Ml- (q) G 


3x3 


denotes the (i,j)-th block of the 


inverse of the M(q) given at ( 11 ), and Qj{q,q) ^ 
is defined at Since the position of the quadrotor is 
X = its acceleration can be written as 


denotes the total mass of the links and the quadrotor. An 
intrinsic formulation of the linearized equations on has 
been developed in [11]. According to it, the variations from 
the equilibrium can be written as 


j u 

=- B{q)u, (30) 

where T{q,q) G and B{q) G 
Assumption 1: The matrix B{q) is invertible for any con¬ 
figuration qi G of the tether chosen such that x = 

YUi=i hQi ^ Dx. 

This assumption is justified by the fact that there is no 
restriction on the acceleration of the quadrotor in and 
therefore it can be arbitrarily changed by the control force 
u, according to Newton’s second law of motion. When the 
quadrotor is on the boundary of i.e., when all of qi is 
identical such that the tether is taut, the control input cannot 
generate any acceleration along qi, due to the constraints 
that the total length of the tether is fixed. This can also 
be observed from the fact that when all of qi are identical, 
the matrix B{q) has a null space spanned by qi, i.e., the 
component of the control force u parallel to qi does not 
affects the acceleration x when all of qi are identical. 
Define a desired trajectory yd{t) as 

yd{t) = x{0)e~'>^ + (1 - <^)(1 - e~'^^)xd, (31) 

for 7 > 0. This satisfies yd{t) = x(0) and limt^oo = (1 “ 
6)xd. Also, yd{t) G Dx for any t > 0, if x(0) G Dx due to 
the convexity of Dx. In other words, yd{t) corresponds to a 
parameterized line connecting the initial point x(0) and the 
intermediate point (1 — 6)Xd. 

Let the position tracking error be = x — yd G 
The control input is designed according to output feedback 
linearization as 

— B {q) q) kx6-x kx^x yd\ ? (32) 

for positive gains kx^kx. 

Proposition 3: Consider the simplified tethered quadrotor 
model described by ([^, where n> 1 and u can be arbitrarily 
selected, with Assumption 1. If the initial position a;( 0 ) 
belongs to the domain given at ( |28l >, then there exist 
controller gains kx^kx such that x{t) G Dx for all t > 0, 
and {cx^Cx) = (0,0) is exponentially stable. 

Proof: See Appendix D. ■ 

2) Stabilization for Tether: The above tracking control 
system guarantees that the quadrotor is translated into the 
intermediate point {1 — 5)xd that is arbitrarily close to the 
actual desired point Xd- But, it does not guarantee that the 
motion of the tether is asymptotically damped out. Therefore, 
we introduce another control system that stabilizes the tether 
as well as the quadrotor. Due to the high degrees of under¬ 
actuation, it is designed based on the linearized dynamics. 

At the desired equilibrium configuration, we have qi = 
— 63 , uji = 0 and Ud = —rriTge^, where ttit = 


X = - y] Mtj{q)ligj{q,q) - 






=-exp(e|i)e3, uj^ = eSoJi, (33) 

where Suii € with • 63 = 0 and 6ui • 63 = 0. This 
yields the following infinitesimal variation 6qi = —^i x 63 . 

Substituting this into ( p3] ) and ignoring the higher 
order terms, the linearized equations can be written as 


Mx + Gx = BSu, (34) 


where x = • • •; ^ corresponds 

to the state vector of the linearized dynamics with C = 
[ei, 62 ] G ei = [1,0,0]^, and 62 = [0,1,0] G The 

control input for the linearized dynamics is Su = u — Ud ^ 
and the matrices M, G G M 2 nx 2 n^ g ^ ]^2nx3 
defined as (see Appendix E) 


Miil2x2 M12I2X2 ••• Minl2x2 

M21I2X2 M22I2X2 ••• ^ 2 n^ 2 x 2 


Mnll2x2 Mn2l2x2 ’’’ ^nn^2x2 


G — diag[(mT — Mgjglil2x2, • • •, (^t — Mg^)glnl2x2], 
~-hC^ei 
—hC^es 


-InC^h 


The control input is designed as 


u = -KxX - KxX - rriTges, (35) 

where the controller gains, Kx,Kx G selected 

such that the linearized dynamics ( [34| ) becomes Hurwitz. 
This provides asymptotic stability of the desired equilibrium 
according to the Lyapunov indirect method. 

In short, the tracking control for the quadrotor position, 
( [^ is engaged first such that the quadrotor becomes suf¬ 
ficiently close to the desired equilibrium. Then, the linear 
control ( [ 35 ] ) is applied to asymptotically stabilize both the 
quadrotor and the tether. 


B. Full Dynamic Model (n > 1) 

The procedures to extend the above control systems into 
the full dynamic model, that incorporate the attitude dynam¬ 
ics of the quadrotor, are identical to described at 

Section III-C. 


C. Numerical Example 

The properties of the quadrotor and the tether are identical 
to Section IIII-DI The initial conditions are chosen such that 
the quadrotor is at x = [2.46, 0, —0.43]^ G and the tether 
is hanging while minimizing the gravitational potential. The 
intermediate position is chosen with 5 = 0.01, 7 = 1, 
and the switching occurs at f = 3. The corresponding 
numerical results are illustrated at Figure where it is 
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Fig. 5. Numerical results for n = 5 (animation available at http:// 
fdel.seas.gwu.edu/CDC15_Fig5.movl 


shown that the quadrotor is translated to the desired position 
asymptotically. In contrast to Figures 3.(b) and 4.(b), the 
vibration of the tether is effectively eliminated at Figure 5.(d) 
and the presented animation. 

In short, the control system presented in this section is 
developed for the flexible cable model when n > 1 at the 
cost of increased complexity. As illustrated by numerical 
examples, the undesired lateral vibrations of the tether, 
observed at Section III are eliminated. The development of 
the dynamic model and the control system design for the 
tethered quadrotor with flexible tether have ben unprece¬ 
dented. While the proposed approach is developed for a 
stabilization problem where the desired link direction is 
flxed, but it is readily generalized to the tracking problems. 

Appendix 

A. Euler-Lagrange equations 

Here, we develop the Euler-Lagrange equations for the 
Lagrangian given by ( p^ and ([^. The Lagrangian is inde¬ 
pendent of R. The derivatives of the Lagrangian with respect 
to O, O are given by 

n 

= Tiq.C = Mgjigez, T>nC = JQ,. 

i=i 

(36) 

Substituting 6R = Rfj into the attitude kinematic equations 
© and rearranging, the variation of the angular velocity can 
be written sls 6Q = r]-\-Qxr] [6]. For the variation model of g^ 
given at we have Sqi = ^iXQi and ii = 6 [7]. 

Let 0 = C dt be the action integral. Using the above 
expression for the variations, and integrating by parts, the 


variation of the action integral can be written as 
<5© = ^ ' {-^Dn£ -nx Dncj ■ V 

+ ^ ^ I ~qi X X dt. 

i=l ^ ^ 

The total thrust of the quadrotor with respect to the inertial 
frame is given by u = —fRe^ G and the total moment 
of the quadrotor is M G with respect to the body-flxed 
frame. The corresponding virtual work can be written as 


ntf 

dW = / u • X g^) + M • 77dt. 

’^^0 i=l 

According to the Lagrange-d’Alembert principle, we have 
6(5 = —6W to obtain 

-f U X = M, 

dt 


-Qi X + gi X J:>q.jC = -kqiU. 


The first equation yields Multiplying both sides of the 
second equation with qi, and substituting ([3^, (36), 


-Qi L]+ OiMg.kgez = -hq^u. (37) 
i=i 

Since qi ' qi = 0, it follows that gi • g’i + ||gi|P = 0. Thus, 
-QiQi = ihx3 - qiqj)qi = qi + 11 * 1 /®- 


Using this identity, © can be rewritten as 

n 

Muqi - qi ^ijQj + MiiWqifqi + qfMgJigez = -q^ku 
i=i 

which corresponds to (0. Rearranging © with the fact that 
iii = -qiCoi - and qiqi = -qfcoi = Coi [7], we obtain 

•GD- 


B. Proof of Proposition 1 

Define an error function, T^(g) = 1 — g • g^. For a positive 
constant f)q < 2, define the following open domain contain¬ 
ing the zero equilibrium, Dq = {(g, cc) G x | < 7 /;^. 

Then, it is shown that 


2 


<^a < - 

- ^-2-^q 


where the upper bound is satisfied for any q e Dq [9]. Define 
a Lyapunov function as 


v, = 


+ c„e, 


q^uj 




which is bounded as 


^q ^q^q — '^q — ^q Pq^qi 

where 2;^ = [||eg||, ||e^||] G M^, and the symmetric matrices 
Pg,Pg G are defined as 


1 

1 

0 , 
1 

CN 

, Pq = - 

2/Cq 

Cq 

2 

_-Cq 1 

’ ® 2 

Cq 1 





































The time-derivative of Vq along (22) can be written as 

Vq = -ZqWqZq, 
x2 • 


where the matrix Wq G 




is defined as ([9]) 


Cq kq 

Cq k(jij 


Cq 


If the constant Cq is sufficiently small, all of the matrices 
Pq^Pq^ Wq are positive-definite, which follows that the zero 
equilibrium of the tracking errors is exponentially stable. 
Substituting ( [ 20 | ) into ( p^ yields T = Td. 


C. Proof of Proposition 2 

The proof is based on the singular perturbation theory [12], 
i.e., if the attitude dynamics is sufficiently fast, the stability 
properties of the reduced system summarized by Proposition 
1 holds. More explicitly, the boundary-layer system corre¬ 
sponds to the attitude dynamics of the quadrotor, and the 
attitude tracking errors exponentially converge to zero at 
the rate proportional to ^ [10, Proposition 2]. The reduced 
system represents the dynamics of the link when R = Rc, 
and from ( |^ , ( [^ , and ( |^ , 

u u 

-fRez = {u • Rcez)Rcez = {u ■ t]—[ 7 ) 7 ]—f? = u. 

I|w|| ll^ll 

Therefore, the reduced system corresponds to the simplified 
dynamic model analyzed at Proposition 1. Then, according to 
Tikhonov’s theorem [12, Thm 9.3], there exists e* > 0 such 
that for all e < e*, the origin of the full dynamics model is 
exponentially stable. 


D. Proof of Proposition 3 

We first show that there exist controller gains such that 
x{f) G Dx for all t > 0. Substituting (32) into (30), we 
obtain the following linear error dynamics. 


Cq* - JvxCx biqr-Cx 


(38) 


Define 




Its time-derivative along the solutions of ( [38] ) is given by 
^ which follows 


kx II 


Therefore, ||ea; 


<U{t)<U{0) = l\\U0)f 


< 


— ^2^ II 

||x||-11^,11 <||x-^,11 = ||e.||. 


for any t > 0. Since 


ikwii<M)ii + 




II4(0)||. 


As Vdit) lies in always, || 2 /d(i)|| < Ei=i 

therefore, if kx is sufficiently small, the above inequality 
guarantees Mm < h as well. Therefore, x{t) G 
for all t > 0 . 

According to Assumption 1, the control input ^32) is well- 
defined, and it is straightforward to show the exponential 
stability of the linear error dynamics given by ( [3^ . 


E. Linearization 

The perturbation model given at ( [^ yields dqi = —^i x 
63 . Substituting it into ([^, dqi is given by 

= iiX -es = X -63 + 0 X X -63) = 6 uji x -63. 

Since both sides of the above equation is perpendicular to 
63, this is equivalent to 63 x x 63) = 63 x { 6 uJi x 63), 
which yields 

^ - (es • ^)e3 = 5 uJi - (03 • SuJi)e 3 . 


Since Ci * ^3 = 0^ we have ^ • 03 = 0. Also, 03 • 6uji = 0 
from the constraint. Substituting these to above, we obtain 
the linearized equation for the kinematics equation: 

= (39) 

Substituting these into ( [T^ , and ignoring the higher order 
terms. 


n 

MiiduJi — MijC^Sujj + Mg.lige\^i = —liC^du + kruTge^^i. 
i=i 


Let C = [ 01 , 02 ] ^ Multiplying the both side of the 

above equation by and rearranging it with the facts that 

( 7^03 = —C^, 03 = diag[— 1 , — 1 , 0 ], C^e^C = —I 2 and 
03 CC^ = 03 , we obtain ([^. 
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